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<^ ■ Abstract 

I investigate the caustics produced by the fall of collisionless dark matter 
in and out of a galaxy in the limit of negligible velocity dispersion. The 
outer caustics are spherical shells enveloping the galaxy. The inner caustics 
are rings. These are located near where the particles with the most angular 
cn ■ momentum are at their distance of closest approach to the galactic center. 

The surface of a caustic ring is a closed tube whose cross-section is a D_4 



catastrophe. It has three cusps amongst which exists a discrete Z3 symmetry. 
A detailed analysis is given in the limit where the flow of particles is axially 



o 

(N 

(N 

^^ ' and reflection symmetric and where the transverse dimensions of the ring are 

^^ ■ small compared to the ring radius. Five parameters describe the caustic in 

■^ ■ that limit. The relations between these parameters and the initial velocity 

Q,. distribution of the particles are derived. The structure of the caustic ring is 

JL \ used to predict the shape of the bump produced in a galactic rotation curve 

H ! by a caustic ring lying in the galactic plane. 
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I. INTRODUCTION 

There are compelling reasons to believe that the dark matter of the universe is consti- 
tuted, at least in part, by non-baryonic coUisionless particles with low primordial velocity 
dispersion [|l| . Such particles are called cold dark matter. The leading candidates are axions 
and weakly interacting massive particles (WIMPs). Before the onset of galaxy formation 
but after the time tgg of equality between matter and radiation, the velocity dispersion of 

the cold dark matter candidates is very small, of order 5va{t) ~ 3 ■ 10^^"^ ( ^^m'^^ ) (^) 

for axions and 6vw{t) ~ 10^^^ (^^) (^) ^'^^ WIMPs, where to is the present age 
of the universe and rria and rriw are respectively the masses of the axion and the WIMP. 
The above estimates of the primordial velocity dispersions, 6va and Sw, are very crude but 
our point is only that in the context of this paper, and of galaxy formation in general, Sva 
and 6vw are entirely negligible. Massive neutrinos, on the other hand, have primordial ve- 
locity dispersion Sv„{t) ^ 5.3 10~^ (^) (^) which is comparable to the virial velocity 
in galaxies and is therefore non-negligible in the context of galaxy formation . For this 
reason, massive neutrinos are called "hot dark matter" . 

Before the onset of galaxy formation, the coUisionless dark matter particles lie on a thin 
3-dimensional (3D) sheet in 6D phase-space. The thickness of this sheet is the primordial 
velocity dispersion Sv. If each of the aforementioned species of coUisionless particles is 
present, the phase-space sheet has three layers, a very thin layer of axions, a medium layer 
of WIMPs and a thick layer of neutrinos. The phase-space sheet is located on the 3D 
hypersurface of points {r,v) : v = H(t)f+ Av{f,t) where H(t) = ^ is the Hubble 
expansion rate and Av{f, t) is the peculiar velocity field. Fig. 1 shows a 2D section of 6D 
phase-space along the (2;, i) plane. The wiggly line is the intersection of the 3D sheet on 
which the particles lie in phase-space with the plane of the figure. The thickness of the line 
is the velocity dispersion 5v, whereas the amplitude of the wiggles in the line is the peculiar 
velocity Av. If there were no peculiar velocities, the line would be straight since z = H{t)z 
in that case. 

The peculiar velocities are associated with density perturbations and grow by gravi- 
tational instability as Av ~ t^'^. On the other hand the primordial velocity dispersion 
decreases on average as Sv ~ t^^/'^, consistently with Liouville's theorem. When a large 
overdensity enters the non-linear regime, the particles in the vicinity of the overdensity fall 
back onto it. This implies that the phase-space sheet "winds up" in clockwise fashion wher- 
ever an overdensity grows in the non-linear regime. One such overdensity is shown in Fig. 1. 
Before density perturbations enter the non-linear regime, there is only one value of velocity, 
i.e. one single flow, at a typical location in physical space, because the phase-space sheet 
covers physical space only once. On the other hand, inside an overdensity in the non-linear 
regime, the phase-space sheet covers physical space multiple times implying that there are 
several, but always an odd number of, flows at such locations in physical space. 

At the boundary between two regions one of which has n flows and the other n + 2 flows, 
the physical space density is very large because the phase-space sheet has a fold there. At 
the fold, the phase-space sheet is tangent to velocity space and hence, in the limit of zero 
velocity dispersion {6v = 0), the density diverges since it is the integral of the phase-space 
density over velocity space. The structure associated with such a phase-space fold is called a 



'caustic'. It is easy to show (see section II) that, in the limit of zero velocity dispersion, the 
density diverges as (i ~ ^ when the caustic is approached from the side with n + 2 flows, 
where o is the distance to the caustic. If the velocity dispersion is small but non-zero, the 
divergence is cut off so that the density is no longer infinite at the caustic but merely very 
large. 

In discussing this type of phenomenon, it is useful to adopt a parametrization of the 
phase-space sheet in the limit of zero velocity dispersion, by giving each particle a 3- 
parameter label a = (01,02,03). The phase-space sheet location at time t is specified 
by the map a — >■ x{a, t) where x is the position in physical space of particle a at time t. 
The velocity of particle a is t7 = % (a, t). The density diverges wherever the map a — > a; 

is singular, i.e. where the Jacobian D = det (||j vanishes. Thus caustics are associated 
with zeros of D. Since D = is one condition on three parameters, caustics are generically 
two dimensional surfaces. 

Y. Zel'dovich emphasized the importance of caustics in large scale structure formation, 
suggesting the name "pancakes" for them. The reason why galaxies tend to lie on surfaces 
1^, such as "the Great Wall", is undoubtedly that the 3D sheet on which the dark matter 
particles and baryons lie in phase-space acquires folds on very large scales, producing caustics 
appropriately called "Zel'dovich pancakes" . 

S. Tremaine recently used the techniques of Catastrophe Theory to catalogue the 
caustics which may occur in observations of structure formation. 

We saw above that where a localized overdensity is growing in the non-linear regime, 
the line which is at the intersection of the phase-space sheet with the [z, z) plane winds 
up in clockwise fashion. The onset of this process is illustrated in Fig. 1. Of course, the 
picture is qualitatively the same in the {x, x) and {y, y) planes. In this view, the process of 
galactic halo formation is the winding up of the phase-space sheet of coUisionless dark matter 
particles. When the galactic center is approached from any direction, the local number of 
flows increases. First, there is one flow, then three flows, then five, seven... The number of 
flows at our location in the Milky Way galaxy has been estimated to be of order 100 |0. The 
boundary between the region with one (three, five, ...) and the region with three (five, seven, 
...) flows is the location of a caustic which is topologically a sphere surrounding the galaxy. 
When these caustic spheres are approached from the inside the density diverges as c? ~ -^ 
in the zero velocity dispersion limit. I call these spheres "outer" caustics to distinguish them 
from the "inner" caustics which are the main topic of this paper. 

To see inner caustics, let us first discuss the case where the overdensity is spherically 
symmetric and all the dark matter particles carry zero angular momentum with respect 
to the center. All particles move on radial orbits then. As a result the galactic center is a 
caustic point where the density associated with each flow in and out of the galaxy diverges as 
c? ~ ^ in the limit of zero velocity dispersion, where r is the radial coordinate. We remarked 
earlier that a caustic is generically a surface. Thus we find a caustic point only because we 
are assuming spherical symmetry and purely radial orbits. If these assumptions are relaxed, 
the caustic at the galactic center will spread over some size a. The notion of caustic point 
is nonetheless useful provided a is small enough. Indeed, for r » a the density will still 
behave as c? ~ ^ and the phase-space structure will be qualitatively the same as for a = 0. 

The question arises what is the inner caustic in the absence of spherical symmetry and 
in the presence of angular momentum. I. Tkachev, Y. Wang and I found P,^ that the 



inner caustic is a ring, i.e. a closed line. The density diverges as (i ~ - in the zero velocity 
dispersion limit where a is the distance to the line |0. However, as remarked earlier, a 
caustic is generically a surface. So a caustic line must again be a special degenerate case. 
Indeed, we will see below that a caustic ring is more precisely a closed tube with a special 
structure. Nonetheless, when the transverse dimensions, called p and q below, are small 
compared to the radius of curvature a of the tube, (i ~ - for a » a » p,q. Hence, it 
makes sense to think of the caustic ring first as a closed line and then, on closer inspection, 
as a closed tube. 

The next question is what is the structure of the tube. We find below that its transverse 
cross-section is a closed line with three cusps, one of which points away from the galactic 
center (see Figs. 5 and 6). In the language of Catastrophe Theory such a singularity is 
called a D_4 catastrophe |p. It has a triality (i.e. a Z^ invariance) which is reminiscent of 
the triality of the Lie group D4, also called S0(8). 

The existence of caustic rings of dark matter results from only two assumptions (see 
section III for a derivation): 

1. the existence of collisionless dark matter 

2. that the velocity dispersion of the infalling dark matter is much less, by a factor ten 
say, than the rotation velocity of the galaxy. 

Only the second assumption requires elaboration. We noted earlier that velocity dispersion 
smoothes out caustics. The question is when is the velocity dispersion so large as to smooth 
caustic rings over distance scales of order the ring radius a, thus making the notion of 
caustic ring meaningless. In ref. this critical velocity dispersion was estimated to be 
30 km/s = 10~^ for the caustic rings in our own galaxy, whose rotation velocity is 220 
km/s. 10~^ is much less than the primordial velocity dispersion 6v of the cold dark matter 
candidates. However the velocity dispersion Av associated with density perturbations also 
smoothes caustics in coarse grained observations. From the point of view of an observer 
with infinite resolution, the effect of Av on a caustic surface is to make it bumpy. However, 
to an observer with poor spatial resolution, this is the same as smoothing the surface. So 
the question is whether the velocity dispersion Av of cold dark matter particles associated 
with density perturbations falling onto our galaxy is less than 30 km/s. The answer is yes 
with very high probability since the infalling dark matter particles are not associated with 
any observed inhomogeneities, and the velocity dispersion of an inhomogeneity as large as 
the Large Magellanic Cloud is still only 10 km/s. The only way the velocity dispersion of 
the infalling dark matter particles could be as large as say 20 km/s is for these particles to 
be part of clumps whose mass/size ratio is 4 times larger than that of the Large Magellanic 
Cloud. But if that were the case, why did these clumps fail to become luminous? 

I have interpreted the appearance of bumps in the rotation curves of NGC3198 and of 
our own galaxy as due to caustic rings of dark matter. That interpretation makes use of an 
additional assumption, namely that the infall of dark matter particles is self-similar |TO| , pA| ,p|. 



This additional assumption is not used in this paper. The goal is to clearly distinguish 
those conclusions which follow exclusively from the two assumptions listed above from those 
conclusions which require the extra assumption of self-similarity. We are however motivated 
by the possiblility that caustic rings may be observed as bumps in rotation curves and we 



derive the shape that such a bump should ideally have. Gravitational lensing by caustic 



rings [|12[ is another technique by which these structures may be observed. 

One might ask whether caustic rings can be seen in N-body simulations of galaxy for- 
mation. Zel'dovich pancakes are seen [|l^. However, caustic rings would require far greater 
resolution than presently available, at least in a 3D simulation of our own halo. Indeed, the 
largest ring in our galaxy has been estimated to have radius of order 40 kpc 0. It occurs in 
a flow that extends to the Galaxy's current turnaround radius, of order 2 Mpc. To resolve 
this first ring, the spatial resolution would have to be of order 10 kpc or smaller. Hence a 
minimum of 2 ■ ■7Ygj^-^^(2Mpc)^ ~ 7- 10"^ particles would be required to see the caustic ring 
in a simulation of this one flow. However, the number of flows at 40 kpc in our halo is of 
order 10 |^. So it appears that several times 10* particles are necessary in a 3D simulation 
of our halo. This is a strict minimum because it only addresses the kinematic requirement of 
resolving the halo in phase-space, assuming moreover that the particles are approximately 
uniformly distributed on the phase-space sheets. There is a further dynamical requirement 
that 2-body collisions do not artificially 'fuzz up' the phase-space sheets. Indeed 2-body 
collisions are entirely negligible in the flow of cold dark matter particles such as axions or 
WIMPs. As a result, Liouville's theorem is strictly obeyed. On the other hand, 2-body 
collisions are present in N-body simulations and thus Liouville's theorem is violated. As a 
result the velocity dispersion is artiflcially increased in the simulations. This may occur to 
such an extent that the caustics are washed away even if several 10* particles are used. The 
best bet to see a caustic ring would obviously be in a 2D simulation of an axially symmetric 
flow. 

This paper is organized as follows. In section II, we give a general discussion of caustic 
surfaces and caustic lines. Caustic lines are a degenerate case whereas caustic surfaces are 
generic. We distinguish two types of caustic line which we call attached and isolated. The 
former are attached to caustic surfaces whereas the latter are not. Caustic rings are of the 
isolated type. In section HI, we discuss the caustics associated with the infall of collisionless 
dark matter particles onto a galaxy. We focus most of our attention on the inner caustic 
rings. We show that the presence of these rings follows only from the two assumptions listed 
above. We flnd that the ring is a tube. Inside the tube are four flows whereas outside the 
tube are two flows. (In addition there is an odd number of flows which are not associated 
with the caustic ring.) In section IV, we give a detailed analysis of the caustic ring under 
the additional assumptions that the flow is axially and reflection symmetric and that the 
transverse dimensions of the ring, p and q, are small compared to the ring radius a. In 
section IVA, we show that under these assumptions the flow near the ring is described in 
terms of flve parameters: a, b, tq, u and s. In section IVB, we relate these parameters to 
the velocity distribution of the infalling dark matter particles. In section IVC, we give a 
qualitative description of the flow of particles on distance scales of order the ring radius a. 
In section V, the results of section IV are used to derive the shape of the bump a caustic 
ring causes in a galactic rotation curve if the ring lies in the galactic plane and the rotation 
curve is measured in the galactic plane. Section VI summarizes the conclusions. 



II. CAUSTICS IN GENERAL 

Consider a flow of collisionless particles with zero (or negligible) velocity dispersion. 
Being collisionless, the particles obey Liouville's theorem. Since they have negligible velocity 
dispersion, the particles lie on a time-dependent 3D sheet in 6D phase-space. The flow is 
completely specified by giving the spatial coordinates x{d, t) of the particle labeled a at 
time t, for all a and t. The 3-parameter label a is chosen arbitrarily. Let dj{x,t), with 
j = l...n, be the solutions of a; = x{a, t). n is the number of distinct fiows at x and t. The 
total number of particles is: 



J daida2da3 J p{ daida2da3 | det f ||] U,{x,t) 

The density of particles in physical space is thus: 

" d^N 1 

^(^'^) = E , , , («.(^,i)) TTu^^i ' (2.2) 

^ daida2da3 \ D{a,t) \sj{x,t) 



where 




D{a,t) = det[-^] . (2.3) 



The formula for the density is, of course, reparametrization invariant. Caustics occur wher- 
ever D = 0, i.e. where the map a ^ x is singular. At the caustic the density diverges. The 
divergence is cut off if the velocity dispersion is finite. 

Generically the zeros oi D are simple, i.e. the matrix || has a single vanishing eigenvalue. 
The condition that one eigenvalue vanishes imposes one constraint on the three parameters 
a. Hence a caustic is generically a 2D surface in physical space. 



A. Generic surface caustics 

Consider a generic surface caustic at a fixed time t. We may reparametrize the fiow 
near the caustic, a ^ P = P{d,t), such that the caustic surface is at P^ = 0. Also, 
in a neighborhood of a point on the surface, choose Cartesian coordinates such that z is 
perpendicular to the surface whereas x and y are parallel. We have then: 

near that point. The 2-dim. matrix Q,fg ^ is non-singular. Since D = at P^ = 0, we have 

z = zo + BPl (2.5) 

for small P^. We may orient the i-axis in such a way that B > 0. Then 



D = 2^B{z-zo) det (^^^) for z > Zo . (2.6) 



Hence, near a caustic surface located at z = zq, the density diverges as -7== on one side of 
the surface. 

Fig. 2a shows a 2D cut of phase-space along the {z, z) plane. The particles lie on a 
line which is at the intersection of the phase-space sheet with the (z, i) plane. The label 
/Ss gives the position of the particles along the line. In the example of the figure, there are 
two caustic surfaces, one dX z = zi and the other ai z = Z2- The two dimensions {x and y) 
into which the caustic extends as a surface are not shown. The density d{z), shown in Fig. 
2b, diverges as —r^ — for z ^> z^ with z > z^ and as —r^ — for z ^ zo with z < Zo- For 
Zi < z < Z2 there are three flows (n = 3) whereas for z < Zi and z > Z2 there is only one 
flow (n = 1). The phase-space sheet "folds back" at z = Zi and z = Z2. That is why the 
map f3 ^ z is singular at these locations. 



B. Line caustics 

Next, consider places where D has a double zero, i.e. where || has vanishing eigenvalues 
for two different eigenvectors. The condition that two eigenvalues vanish defines a line. This 
line is the location of a more singular kind of caustic which we call 'line caustic'. 

In a small neighborhood of a line caustic let us reparametrize the flow {a — > /3) such 
that /9i = /?2 = defines the line in the new coordinates. Also, choose Cartesian coordinates 
such that z is parallel to the line and x and y are perpendicular to it. In this neighborhood 
we have: 

with -M- 7^ 0. At /9i = /?2 = 0, the matrix g/gg -. vanishes since it is 2x2 and has vanishing 
eigenvalues for two different eigenvectors. Thus for small f3i and f32, 

x = xo + ^XnPl + Xuf3if32 + \x220i + 0(/9') 

1/ = 1/0 + ]^xxfil + Yx2^xfii + ]^22fil + ^{^f) • (2.8) 

Linear combinations x — xq + 7(1/ — yo) are complete squares provided: 

l'{YxxY22 - nl) + 7(>^ii^22 + Xnr22 - 2X12F12) + XnX22 - Xl^ = , (2.9) 

which has two solutions: 

7± = ^,^ ^ 772t{ —^11-^22 — ^^11^22 + 2Xi2yi2 

± [(1^11X22 + XnF22 - 2X12F12)' - 4 (rnl^22 - Y,\) (XnX22 - Xf^)] '^'} . (2.10) 

If the denominator ^11^22 — ^12 — 0; '^^^ should interchange the role of x and y. If Y11Y22 — 
Y^2 = and X11X22 — XI2 = 0, a; — xo and y — yo are complete squares to begin with. Note 
that Im7-|- may differ from zero and that 7_ = 7I in this case. We have: 



1 2 

x± = X - xo + -f±{y - yo) = -X± {pi + r]±p2) (2.11) 

with 

X± = Xn+7±>'ii (2.12) 

and 

-^12 + 7±^12 /^ -ION 

^± = y , 7F- ■ (2-13) 

Xii + 7±yii 
In terms of these new quantities: 

^-^det(^^]=2./5^;]^:^.^-±^^ . (2.14) 

\d{Pup2)J ^ 7+-7- 

There are two cases to consider, depending on whether Ini7-|- 7^ or Im7-|- = 0. 

If Im7-|- 7^ 0, let 7-1- = 7i ± ^72, X± = Xi ± 2X2, x± = Xi ± 1x2 and r^-t = r^i ± iri2. 
7i,72,Xi,X2,Xi,a;2,?7i and rj2 are real. We have: 



V2 
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D2 = 2 ^ ^Xl + X| (x - xo + 71(1/ - 1/0))' + 72(1/ - Vof ^'" • (2.15) 

72 ^ L J 

1^2 7^ everywhere except at (x, y) = (a;o,t/o). The density c? ~ - with 

o- = \l{x - xo + 71(1/ - 1/0))^ + 72(1/ - Vq? ■ (2.16) 

An illustrative example of a Im7-|- 7^ line caustic is: 

x = x^ + x]^{Pl-Pl) , y = t/o + F/9i/92 (2.17) 



for which D = 2JY'^(x — xq)"^ + X^(y — t/o)^- We call a Im7-|- 7^ line caustic "isolated" 
because it is not attached to any surface caustic. Ring caustics associated with the infall of 
collisionless dark matter onto a galaxy are line caustics of the isolated type. They will be 
discussed in detail in sections III and IV. 

If Im7-|- = 0, then x±, X± and ri± are real. From Eq. (2.11) and Eq. (2.14) we see that 
D has simple zeros along two surfaces, x+ = and x_ = 0, and a double zero on the line 
where these two surfaces meet. Thus the Im7-|- = line caustic is "attached" to two generic 
surface caustics. The following is an illustrative example: 

x = xo + ^Xf3l y = y, + lYl3l. (2.18) 



In this case D = 2JXY{x — XQ)(y — yo), and hence the density d ~ (^ ^^oj (y yo) ^^^ X,Y > 

V ^(x-xo){y-yo) 

0. We have no use for the "attached" type of hne caustic in this paper. 



III. CAUSTIC RINGS 

In this section we give a general discussion of the caustic rings that are associated with 
the infall of colhsionless dark matter onto a galaxy. No symmetry is assumed. We do assume 
that the velocity dispersion of the dark matter particles is small compared to the rotation 
velocity of the galaxy. We set the velocity dispersion equal to zero and derive the existence 
of caustics. A small velocity dispersion provides a cutoff, so that the density at the caustic 
does not become infinite but merely very large. 

Fig. 3 shows successive time frames of a set of colhsionless particles falling through a 
galaxy. The particles move purely under the effect of gravity. In Figure 3a, the particles 
are at first 'turnaround', i.e. they are about to fall onto the galaxy for the first time in 
their history. They are located on a closed 2D surface surrounding the galaxy in physical 
space, called the "turnaround sphere" . Fig. 3 shows the intersection of this sphere with the 
plane of the figure, as it evolves in time. For the sake of definiteness it is assumed that the 
particles carry net angular momentum about the vertical axis. Qualitatively speaking, the 
turnaround sphere is spinning about the vertical. The particles near the top (bottom) of 
the sphere in frame a carry little angular momentum and end up near the bottom (top) of 
the sphere in frame f after falling through the galaxy. The particles near the equator carry 
the most angular momentum. They form a ring whose radius decreases in time down to 
some minimum value, reached near frame d, and then increases again. Generally the radius 
of the sphere at second turnaround is smaller than at first turnaround because the galaxy 
has grown by infall in the meantime. After second turnaround the sphere falls back in and 
repeats the same qualitative sequence till third turnaround, and so on. 

There is a generic surface caustic associated with the mth turnaround where m = 
2,3,4,.... These caustics are located near where the particles of a given outfiow reach 
their maximum radius before falling back in. To see this, parametrize the fiow at a given 
time t by x{a, P, to; t) where to is the time the particle was at last turnaround and a, P - for 
example, spherical coordinates - tell us where the particle was on the turnaround sphere at 
that time. We define: f° = ^, a;^ = |^, a;^ = ||, and x = ^. In terms of these: 

D= fO ■ (x^xx^) . (3.1) 

Since the discussion is for a fixed time t, let us not show the t dependence explicitly further. 
Let us assume that x{a,(3,tQ) is near particles at their mth turnaround with m = 2,3,4... 
The spheres {x{a, /3, tg) : Va, j3} with tg considerably larger or considerably smaller than to 
are inside the sphere {x{a,(3,to) : \/a,P}. This implies that for all a,P there is a to(a,/5) 
such that x^{a, P,to{a, P)) is parallel to the sphere {x{a', P',to{a, P) : Wa',P'}, and hence 
where D = 0. Therefore the sphere {x{a, P,to{a, P)) : \/a,P} is the location of a generic 
surface caustic. These outer caustics were mentioned in the introduction as the boundaries 
between the regions with n fiows and those with n + 2 fiows where n = 1,3,5... We do 
not discuss them further in this paper. Our focus is upon the inner caustics which, as we 
will soon see, have the shape of rings (closed tubes) and which are located near where the 
particles with the most angular momentum in a given infiow reach their distance of closest 
approach to the galactic center before moving back out of the galaxy. 

To see inner caustic rings, let us return to the sphere of Fig. 3. During each infall-outfall 
sequence, the sphere turns itself inside out. Indeed, a particle which is part of the fiow 



and is just inside the sphere in frame b of Fig. 3 is outside the sphere in frame e, and 
vice- versa. There is therefore a ring of points in space-time which are inside the sphere last. 
The intersection of this space-time ring with the plane of the figure is at two space-time 
points, one located at the cusp in frame d, the other at the cusp in frame e. Since there is 
a continuous flow of spheres falling in and out, the ring just defined is a persistent feature 
in space. We now show that it is the location of a caustic. 

Fig. 4 shows the infall sphere in a neighborhood of space and time where it is completing 
the process of turning itself inside out. We choose Cartesian coordinates such that y is 
parallel to the ring at that point; x and z are as shown. We parametrize the flow in a small 
neighborhood of the ring by cJ = (q;,/3, to) such that f| = |f = 0- As before, to labels 
successive infall spheres and may be taken to be the time of their last turnaround. Thus: 

jj^^ (^ d^ _dx_ dz\ 

dp \da dto dto da) ' ^ ' ^ 

Fig. 4a shows an infall sphere just before it reaches the ring. The line is the intersection 
of the infall sphere with the plane of the figure, a labels points along the line. We have 
1^ = at points A and B because the line is parallel to the x-axis at these points. Similarly, 
1^ = at point C. Fig. 4b shows an infall sphere at the moment it reaches the ring. Points 
A, B and C in Fig. 4a have moved to point E in Fig. 4b. Hence |f = f^ = and therefore 
D = aX point E, which is thus the location of a caustic. 

In general, D only has a simple zero at E, which means that E is the location of a surface 
caustic. We show below that the complete inner caustic surface is a closed tube. E is a 
point on this tube. In the limit where the transverse size p of the tube goes to zero, both 
1^ = 1^ = and ^ = ^ = at E. D then has a double zero at E. In that limit the tube 
caustic becomes an isolated line caustic where the density diverges as -. When p is finite 
but small, the density d{a) ~ - for a » p but has a more complicated behaviour for a ^ p. 

Let us show that the inner caustic must have the topology of a tube, starting with the 
case of Fig. 3 and generalizing from there. In Fig. 3b, aP (not shown explicitly) is everywhere 
pointing outward of the infall sphere because later infall spheres are outside this one. In 
Fig. 3f, of is pointing inside because later infall spheres are inside this one. This implies 
that during the infall of either vanishes at some space-time points or becomes parallel to 
the sphere. Since D = at such points, they are the location of caustics. Since D is a 
continuous function of a, the caustic must lie on a closed surface. Now, if we follow the 
motion of points which are near the top (bottom) of the sphere in Fig. 3b and end up near 
the bottom (top) in Fig. 3e, of always points up (down), a* does not vanish and is not 
parallel to the sphere at any time between these two frames for these points. This implies 
that there are no inner caustics within some cylinder extending from top to bottom in the 
spatial volume under consideration in Fig. 3. On the other hand, for the points near the 
equator in Fig. 3, a* is pointing outward during infall and is pointing inward during outfall. 
Thus if we track a point near the equator, at some time aP either vanishes or is parallel 
to the sphere. The points where this happens lie on a closed surface which is outside the 
previously defined cylinder but wraps around it. That surface must therefore be a tube. 
The tube is located near the equator, where the particles with the most angular momentum 
are at their distance of closest approach. 

Consider now the most general angular momentum distribution on the turnaround 
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sphere. The angular momentum is a continuous two-dimensional vector field on the sphere. 
It is well-known that such a vector field must have (at least) two zeros. If we track the 
particles in the neighborhood of these angular momentum zeros, we find that their aP does 
not vanish and is not parallel to the sphere at any time during the infall-outfall sequence. 
Therefore the inner caustic appears only in the fiow of particles which are some distance 
away from both zeros. Since that set of particles has the topology of a closed ribbon, and 
the previously defined caustic ring (the set of points E which are in the turnaround sphere 
last) goes around this closed ribbon once, the inner caustic must be a closed surface with one 
handle, i.e. a closed tube. Note that the closed tube may self-intersect and may have a com- 
plicated shape. However, if the angular momentum is dominated by a smooth component 
which carries net angular momentum, as was assumed in Fig. 3, then the tube resembles a 
circle. If there is no angular momentum at all, the tube reduces to a point at the galactic 
center. As was mentioned in the Introduction, the galactic center is a point caustic in that 
case, where the density behaves as d{r) ~ ^. 

Fig. 5 shows simultaneous (same t) infall spheres near the caustic corresponding to five 
different initial times: toi > ^02 > ••• > ^05 • The five numbered dots show the positions 
x{a,tQk),k = 1...5, for some fixed a and the five different initial times. Following the dots 
gives one a qualitative picture of the flow in time. The toi sphere is falling in but has not 
yet crossed itself, as in frame b of Fig. 3. The to2 sphere has crossed itself but has not yet 
completed the process of turning itself inside out, as in frame c of Fig. 3. The tos sphere is 
just completing the process of turning itself inside out, as in frame d of Fig. 3. The cusp 
at point E is the location of a caustic for the reason given earlier. Let aE be the value of a 
at E. Thus x{aE,to3) is the position of E. In the example of the flgure, the particles at the 
cusp are moving to the left. Thus at {aE, tos), x is pointing to the left whereas of is pointing 
to the right of the tos infall sphere. For smaller initial times, such as to5) 2* is pointing to 
the left of the infall sphere. Let to4 be the initial time and F be the point where af{aE,tQ) 
crosses the to sphere. In view of Eq. (3.1), F is the location of a caustic as well. Consider 
any point which is far from both E and F. At such a point there are two flows because the 
sphere passes such a point twice, once on the way in and once on the way out. Consider also 
a point located between E and F. At such a point there are four flows because the sphere 
passes by four times: at flxed t, twice for initial time to between to2 and tos, once between 
to3 and to4 and once between to4 and to5. There is therefore a flnite compact region in the 
plane of Fig. 5 inside of which there are four flows and outside of which there are two flows. 
The boundary of this region, shown as a closed dashed line with three cusps, is the location 
of the tube caustic. The reason for the three cusps will be given in the next section. 

In the limit where the tube shrinks to a line, D has a double zero on the line. Indeed 
E and F coincide in that limit. aP is parallel to the ring at E because of is parallel to the 
infall sphere at F which coincides with E and the infall sphere has a cusp at E. Eq. (3.2) 
shows that D has a double zero then which means that the ring is a line caustic. It is of the 
isolated type since it is not attached to any surface caustics. 

In the example of the flgure, the initial time tos of E is after the initial time to4 of F. 
The opposite is equally possible. Indeed the time reverse of the sequence of Fig. 5 is also 
a possible sequence and it has the times of E and of F in the opposite order. On the other 
hand, the density is invariant under time reversal of the flow. This shows that point F is 
always closer to the galactic center than point E. Note also that if the flow is time reversal 
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invariant (i.e. if the outfall sequence is the time reverse of the infall sequence), the tube 
caustic collapses to a line caustic because to4 = tos and hence E and F coincide. 

If the successive spheres all fall in exactly the same way, the caustic ring is stationary. 
In general, however, the trajectories of successive spheres change with time t, albeit slowly. 
As a result the caustic ring moves about. However, whether or not it moves, the caustic 
ring is perfectly sharp in the limit of zero velocity dispersion. 

IV. THE AXIALLY SYMMETRIC CASE 

We assume in this section that the flow of dark matter in and out of the galaxy is 
axially symmetric about z and symmetric under reflection z -^ —z. We also assume that 
the transverse dimensions, called p and q below, of the ring caustic are small compared to 
the ring radius a. As before, we assume that the dark matter particles are collisionless and 
we neglect their velocity dispersion. We use the following parametrization of the flow. Let 
R(to) be the turnaround radius in the z = plane at time to. Then let ^(^o, V^o, ^o; t) be the 
position at time t of the particle that was at the location of polar coordinates {60, ipo) on the 
sphere of radius -R(to) at time to- The density is given by Eqs. (2.2-3) with a = {9o, ipo, to) 
and n = 2(4) outside (inside) the closed tube caustic described in a general fashion in the 
previous section. 

However, since we assume axial symmetry the motion is in effect two-dimensional. Let 
p{a, to; t) and z{a, to; t) be the cylindrical coordinates at time t of the ring of particles which 
start with azimutal angle 60 = ^ — a ai initial time to. One readily shows, repeating for the 
axially symmetric case the steps which led to Eqs. (2.2-3), that the density is given by: 



1 J^ (fN ^ 
d{p,z,t) = - — 2^ -— — (a, to 



where 



2'Kp ^ dadtf) \ D2{a,to) 



(4.1) 

{a,to)={a,to)j 



and (Q;,to)j, with j = l...n, are the solutions of p = p{a,tQ;t) , z = z{a,tQ;t). Under 
a reparametrization of the flow (a, to) -^ [Q;'(Q;,to), tQ(Q;,to)], the determinant transforms 
according to: 

D,{a, to) = D'.ia', Q det (j^) ■ (4-3) 

In particular, 

D2{a,to)=D',{a',Q (4.4) 

for an a-dependent time shift: a' = a, tg = to -|- Ato(Q;). 
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A. Flow at the caustic 

General characteristics of the flow near the caustic are described in section III. Figure 
5 constitutes a summary. In the reflection symmetric case, points E and F have z = and 
a = 0. Let us reparametrize the flow a -^ a, t^ —>■ t = to + Ato(a) such that z{a, r = 0) = 
for all a. We expand p and z in powers of a and r keeping terms up to second order only. 
Since we assume the reflection symmetry: z —>■ — z, p -^ p, a ^ —a, and r ^ r, we have: 

z = +baT (4.5) 

p = PQ-CT + -ut'^ - -sa^ , (4.6) 

where b,po,c,u and s are constants, po is the p-coordinate of point E. We may rewrite 
Eq. (|;6|) as 

p = a+-M(r-ro)^- -sa% (4.7) 

where tq = c/u and ci = Po — \ ^- ct is the p-coordinate of point F. The parameters 6, u 
and s are positive, h is positive because the flow is from top to bottom (bottom to top) for 
particles with a > (a < 0). u is positive because the particles with a = are accelerated 
outward by the angular momentum barrier, s is positive because at r = 0, the particles 
with a 7^ are at p < po. tq can either be positive or negative, a is always smaller than po, 
i.e. point F is always closer to the galactic center than point E, as noted earlier. 
The determinant is 

D2(«, r) = -b[ur{T - To) + sa^] . (4.8) 

It vanishes for 



a = ±\l^r{To - r) (4.9) 

with < r < To if To > 0, and -tq < r < if tq < 0. Substituting Eq. ( ^^9]) into Eqs. (4.5) 
and (4.7), we find a parametric (r = parameter) representation of the surface of the tube 
caustic: 

p = a + -u{t - Tq) (2r - tq) 



Z = ±bJ-T^{To-T) . (4.10) 

Fig. 6 shows a cross-section of the tube. Its dimensions in the p and z directions are 
p = \uTq and q = ^^ -^ p respectively. The caustic has three cusps. It is shown below 



2""U ""^ -i 4 Vili 

that the appearance of the three cusps is not special to the assumed axial and reflection 
symmetries nor to the fact that we expanded only to second order in Eqs. ( |4.5| ) and ( [4.6| ). 
First let us note that the caustic has a Z^ symmetry after an appropriate rescaling. 
Define the rescaled parameters: 
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T = -, A=\[^-a, (4.11) 

To V n To 



and the rescaled and shifted (in the p direction) coordinates: 



Z = lJ^\z, X = ^{p-a)-l. (4.12) 

b y U Tn UTn 4 



'0 "'0 

In terms of these, Eqs. ( [4. 51) and (|4.7| ) become: 



Z = 2AT, X = {T-lf-A'^ -J. (4.13) 



These relations are invariant under the discrete transformation: 

Z' = -^-Z+^X , X' = -^Z-^-X 
2 2' 22 

T' = -It+^A + 3/A , A' = -y^T-lA + V3/4 



(4.14) 



whose cube is the identity. In the X-Z plane, the transformation is a rotation by 120°. It 
transforms the 3 cusps of the caustic into one another. Hereafter, let us call the shape 
of Fig. 6 a "tricusp". In the language of Catastrophe Theory 0, the tricusp is a D_4 
catastrophe. 

The apparent reason for the two cusps which are not on the p-axis in Fig. 6 is that -^ 
and ^ both vanish for r = |ro. It might be thought that the simultaneous vanishing of 
-^ and ^ is an accident peculiar to our assumptions of symmetry and/or our limiting the 
expansion of p and z to terms of second order in a and r. But this is not the case. Indeed, 
consider the equation for the location of a generic caustic in 2 dimensions: 

It defines a(r) such that \p{r) = p{a{r),T), z{t) = z{a{r), r)] is a parametric representation 
of the caustic location. Wherever p(r) has an extremum, z{t) may be expected to have an 
extremum as well since 

dp dp da dp , , . 



and Eq. ( [4.15| ) imply: 



dr da dr dr 



dz dz da dz , , 

dT da dr dr 



if ^ is finite. The cusp at point E may appear to have a different origin. Its apparent 

3 

reason is that z ~ ±r 2 near r = which is at the boundary of the range of r. However, this 
circumstance is an artifact of the parametrization used. Indeed, we saw that the three cusps 
are transformed into one another by a Z^ symmetry. As a corollary, the cusp on the p-axis 
can be given the same parametrization as the other two. In conclusion, the appearance of 
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three cusps in the cross-section of a ring caustic is not an accidental consequence of our 
simphfying assumptions. Without those assumptions, the cross-section does not have the 
exact shape of Fig. 6 but it still has three cusps, at least for small deviations from the 
description given above. 

Inside (outside) the tricusp, there are four (two) flows. For example, if we restrict 
ourselves to the z = plane, then a = or r = 0. If a = 0, then p = a + |m(t — tqY > a. 
If r = 0, then p = a + ^uTq — |sa^ < po- Thus for z = and p > a, there are 2 flows: 

a = , r = ro±J-{p-a) , (4.18) 

V u 

for which 



^ = , 1^ = ±^2u{p - a) . (4.19) 

We may call these the "in and out" flows. For z = and p < Pq, there are 2 other flows: 

r = , a = ±^^(po-p), (4.20) 

for which 



|i = ±(,y'?(p„-rt . | = -«r,. (4,21) 

We may call these the "up and down" flows. Therefore, in the z = plane, there are four 
flows (in, out, up and down) for a < p < po, whereas there are two flows (up and down) for 
p < a, and two flows (in and out) for p > po- 

Away from the cusps, the caustic is a generic surface caustic, as described in section IIA. 
Thus, if one approaches the boundary of the tricusp from the inside and away from any 
of the cusps, the density increases as rf ~ ^ where a is the distance to the boundary. If 
the boundary is approached from the outside, away from any of the cusps, then the density 
remains finite until the boundary is reached. For example, for 2; = and p — a ^ 0+, -D2 — 



=F6roA/2-u(p — a) for the in and out flows, and hence the density associated with these flows 
increases as -7^ — . 

•yp—a 

Near the cusps, the behaviour depends upon the direction of approach. For z = and 
p — Po ^ 0_, D2 — — 26(po — p) for the up and down flows and hence d ~ —^. For z = 

/^2{j,j.2^U|2\ 1/3 

and p — Po -^ 0+, the density remains flnite. For p = po and z ^ 0±, | D2 '~ ■'^ ' ° ' 

for one of the flows and hence d ~ r-kja- 

The tube caustic collapses to a line caustic in the limit Tq ^ with -^ flxed. In this 
limit, p,q = 0, 



us 



D2 = -b{uT^ + sa^) = -2^b^{p-ay + usz^ , (4.22) 

and hence the density 
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dU z) = = . 4.23 

2vrp dadto Jh'^{p - af + usz"^ 



When p and q are finite but much smaller than a, Eq. (4.23) is approximately valid for 
p^q « a « a since the terms of order tq in Eqs. ([4.5|), ( |4.7| ) and (^78|) are small in this 
regime. 

B. Relations between the caustic ring parameters and the initial velocity 
distribution of the dark matter particles 

We saw in section IVA that the flow near the caustic is parametrized by the quantities 
6, a, u, s and tq. In this section we relate these quantities to the velocity distribution of 
the dark matter particles near a = at the initial time to- We assume that the gravitational 
potential of the galaxy is spherically symmetric, so that the angular momentum of each 
particle is conserved. As before, we assume p,q « a. Thus the particles participating in 
the flow at the caustic have | a \« 1. Also the time scale over which the particles cross 
the caustic is short in that limit, which implies that their velocity is nearly constant while 
they do so. 

As before, i?(to) is the turnaround radius of the particles in the equatorial plane (a = 0) 
at time to- Thus the radial velocity of the particles at a = 0, t = to, and r = -R(to), vanishes. 
Consider the velocity distribution of the particles on the sphere of radius -R(to) at time to 
near a = 0: 

v{9o,ipo,to) = (^oviiia) + foVria) - 9ov±{a) 

= ^obii(O) + lvi;{0)a' + 0{a')] + fo[^<(0)a2 + 0{a^)] 

+ 9o[-av'^{0) + 0{a^)] , (4.24) 

assuming axial and reflection symmetry, fo, <^o and 6q are the unit vectors in spherical 
coordinates at location {9o,ipo) on the sphere. The dependence of v\\{a),v±{a) and Vr{a) 
upon to is not shown explicitly but is understood. By definition v'j_{0) = ^(0), ^[[(O) = 

-;^(0), and so on. We need only discuss the motion of the particles initially at y^o = f since 
the motion of particles with differing values of ifQ are trivially related by axial symmetry. 
The particle which is initially at (/Pq = f , 6*0 = | — a has angular momentum: 

/(«) = RfoXv{^- a, |, to) = iia)iia) (4.25) 

where R = R(to) and 



i{a) = R^vf^ia) + vl{a) , (4.26) 

and 

£{a) = cos(l){a){cosa z — sina y) + sin((){a) x , (4.27) 

with 
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(j){a)=tan-^ (^^4^1 ■ (4-28) 

\v\i{a) J 

Let {x' , y') be Cartesian coordinates in the plane of the orbit of that particle, with origin at 
the galactic center. The {x',y',z') coordinates are related to the {x,y,z) coordinates by a 
rotation of angle a about the a;-axis followed by a rotation of angle 0(a) about the y'-axis: 

sina sin(f){a) —cosa sin(f){a) 

cosa sina \ \ V \ ■ (4.29) 

—sina cos(l){a) cosa cos(l){a) 

In the {x',y') coordinates, the particle starts at tq = Ry' with initial velocity VQ{a) = 
— R^' "*" '^r{o!)y' ■ See Fig. 7. At the moment of its closest approach to the galactic center, 
the particle moves with approximately constant velocity: 

x'(Q;,to; i) = —fm cos5 + V {t — tm) sin6 

y'{a, to; t) = -Tm sin5 -V {t- tm) cos6 

z'{a,to;t) = , (4.30) 

where r^ and t^ are the distance and time of closest approach, V is the magnitude of 
velocity then, and 6 is defined in Fig. 7. The quantities r^, tm, V and 6 depend on a and 
to- On the time scale over which the particles cross the ring caustic, the flow is very nearly 
time-independent, i.e. x{a, to; t) = x{a, t — to). Hence, we replace t — t^ by to,™ — ^o in Eqs. 
( [4.30| ) where to,m(«) is the initial time of the particles which are at their closest approach 
at time t. The flow at fixed t is then given as a function of a and to by: 

x'{a,to) = — rm(a) cos6{a) + V{a) sin6{a) (to,m(a;) — tg) 
y'{a,to) = -r.m{<y) sin6{a) - V{a) cos6{a) (to,m(a) - to) 
/(a,to)=0 . (4.31) 

In Eqs. ( [4.31| ) and henceforth, the t-dependence of x',y',rm,S,V and to,m is not explicitly 
shown. 

Using Eqs. ([4.29|) and ( |4.31| ), we have: 



2;(Q;,to) = rjn{(y)[cosa sin(j){a) cos5{a) — sina sin5{a)\ 

+ V{a)(tQ — to,m(a))[cosQ; sin(f){a) sin6{a) + sina cos6{a)] . (4.32) 

As in section IVA, we reparametrize the flow to ^ r = to — ro(a) such that z{a, r) = at 
r = for all a. Thus: 

z{a, t) = V{a) T [cosa sin(j){a) sin6{a) + sina cos6{a)] . (4.33) 

The time shift To (a) is given by: 

V (a) (ro(a;) — to,m(a)) [cosa sin(j){a) sin6{a) + sina cos6{a)] 

+ rm{Q)[cosa sin(j){a) cos6{a) — sina sin6{a)] = . (4.34) 
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Combining Eqs. (^1291) , (WM) and (|;3|), we have: 



x{a, t) = — cos(f) sinS Vt — r^ coscl) 



cosS — sin6 



cosa sincj) cos5 — sina sin6 
cosa sind) sinS + sina cos6 



y{a,T) = — Vr^sina sincp sin6 — cosa cos6) 



sina sincj) cos6 + cosa sinS 



[sina sincj) sin6 — cosa cos6) 



cosa sincj) cos6 — sina sin6 
cosa sincj) sin6 + sina cos6 



(4.35) 



with V = V(a), cj> = cj){a), r^ = ?"m(a) and 5 = 6{a). 

We now compare Eqs. ( |4.33|) and ( [4.35| ) with Eqs. ([4.5|) and ( |4.7| ) to extract b, a, u, tq 
and s. cj>{a) is an odd function of a whereas rm{a), V{a) and 6{a) are even. Thus: 



(a) = cj)'{0)a + - cj)"'{0)a^ + O(a^) 
6 



(4.36) 



r™(«) = r^(0) + - C(0)«2 + 0(«^) , 



and so on. Comparing Eqs. ( [4.5|) and ( [4.33|) , we have 

b = V{0){cos6{0) + cj)'{0)sinS{0)) . 
We may rewrite Eqs. ( [4.35|) as 



a^ + uaT^ = Xo(0) + i/o(0) 

-uaro = xo(0) Xi(0) + yo(0) yi(0) 
ua = xl{0) + yl{0) . 



(4.37) 



(4.38) 



x{a, t) = xo{a) + xi(a)r 

y{a,T) = yo{a)+yi{a)T , (4.39) 

with the appropriate definitions oi Xo{a),Xi{a),yo{a), and yi{a). Therefore 

p{a, rf = x{a, rf + y{a, rf 

= xoiaf + yo{af + 2 {xo{a)xi{a) + yo{a)yi{a)) r + (^Xiiaf + yi^af) r^ , (4.40) 

which is to be compared with the square of Eq. ( [4. 71) : 

p(a, rf = a^ + ua{T - r^f - saa^ + 0(r^, r^ , aV) . (4.41) 

For a = this yields: 



(4.42) 



From Eqs. (|4.35|) and ( [4.39| ), we have 



a;o(0) 
a:i(0) 



-r™(0) 



cos6{0) — sin6{0) 



-V{0) sin6{0) 



yo{0) = -r„(0) 



sin6{0) + cos6{0) 



4>\0) cosS{0) - sinS{0) 
0'(O) sinS{0) + cosS{0) 

0'(O) cos5{0) - sin5{Q) 
(/.'(O) sin5{0) + cos5{0) 



Hence: 



y^(0) = V{0)cosS{0) . 



V{Of 



u 



To 



r„(0) 0'(O) cos5{0) - sin5{0) 
V{0) 0'(O) sinS{0) + cosS{0) 

Comparing Eqs. ([4.40|) and (|4.41| ) for a 7^ 0, we have: 

-sa = {xo{0)x'^{0) + yo{0)y'^m 
A somewhat lengthy calculation yields: 

^^ _ C(o)(i + 0W) 

[(j)'{0)sin6{0) + cos(5(0)]2 

0'(O)cos5(O) - sin5{0) 



+ 



[0'(O)sm5(O) + cos(5(0)]3 

^(0'"(o)-0W)-(i + 0'(of)no) 



I ^'(°) 



Our derivation of the caustic parameters assumes that p and q 
smaller than a. Using Eqs. ([4. 441), we have: 



v^ 



p 



p 



1 U'{0)cos6{0)-sin6{0y 
2^ \(f)'{0)sin6{0) + cos6{0) ^ 



(4.43) 



(4.44) 



(4.45) 



(4.46) 
are much 

(4.47) 



Thus the treatment requires that (f)'{0) and 5(0), or at least the combination (f)'{0)cos6{0) — 
sin6{0), be small compared to one. 

Eqs. ( [4. 381 ), ( [4. 441 ) and ( [4. 461 ) express the caustic parameters in terms of the values of 
(p{a), Tmia), V{a) and 6{a) and their first few derivatives at a = 0. 0(a;) is given in terms 
of the initial velocity distribution by Eq. (|4.28| ). Moreover, angular momentum conservation 
implies: 



rm{a) 






(4.48) 



with i{a) given by Eq. ( [4.26] ). Thus, to achieve our goal of determining the caustic param- 
eters in terms of the initial velocity distribution, it remains to express V{a) and 6{a) in 
terms of ^^(q;), f||(«) and v±{a). 
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This last step can only be carried out if we adopt a model for the galactic gravitational 
potential U{r, t) in which the particles fall. To illustrate the process, let us adopt the time- 
independent potential: 



2 1 i'R\ 



U{r) = <, In ^-j (4.49) 

which yields (perfectly) flat rotation curves with rotation velocity Vrot- Since particle energy 
is conserved for this potential, we have: 



V{a) 



2w^„4ln ( - j + ^^r(a) + vfia) + vl{a) 



(4.50) 



The angle S = 6 {i, v^) was determined numerically for the potential U{r) by solving the 
equations of motion of a particle falling from the initial position tq = Ry' with initial 
velocity vq = ^ x' + Vrjj'] see Fig. 7. Table I gives 5 as a function of j = -^^ and p = -^. 
A simple velocity distribution is that corresponding to an initially rigidly rotating 
turnaround sphere: v^{a) = 0, Vr{a) = and v\\{a) = f|j(0)cosa. Let us explore what 
happens in this case. Since 0(a) = 0, we have: 

p = -a tan^(5(0) 

- ^: , -°'^^(°' . (4.51) 

^ - tan5(l))S"(l)) 

Typical values in fits of the infall model [|n|,^0 to observed properties of our galaxy are 
j ~ 0.25 and a ~ , -^ ' ~ 0.1-R. Table I shows that S is slowly varying and of order 

y/2ln{R/a) 

0.5 then. For 0'(O) = 0, this implies p ~ 0.15 a. For small j,V ~ Vrot\/'^^^{R/o-) is 
approximately a-independent. For the above special velocity distribution, we have then 



fm{oi) = rm(0) cosa and hence q ~ ^^ p cos 5(0) ~ p . 

C. Flow some distance away from the caustic 

In this section, we give a qualitative description of the flow associated with a caustic 
ring on distance scales of order a, the ring radius. To this effect, we choose the special 
initial velocity distribution corresponding to an initially rigidly rotating turnaround sphere: 
Vr{a) = v±{a) = 0, ^11 (a) = v\\{0)cosa. Also we neglect the a-dependence of V and set 
5 = 0. Table I shows that the latter approximation is valid only for j << 1, i.e. a << R. 
We also take the velocity V of each particle to be constant while it travels distances of order 
a. These crude approximations yield a description of the flow which is topologically correct 
and which can be derived by simple analytical methods, but which is likely to agree only 
qualitatively with actual flows. In particular, since (p{a) and 6{a) are set equal to zero, the 
tricusp structure of the caustic is shrunk to a point. 

For the assumed velocity distribution, 

i{a) = irnax cosa (4.52) 
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with Irnax = Rv\\{0). Siiice we neglect the a-dependence of V, 

^m(«) = o cosa (4.53) 

with a = ^^!^. Since we set (j){a) = and 5{a) = 0, Eqs. ( |4.33| ) and ( [4.35| ) become: 



x{a, t) = —a cosa 
y{a, t) = Vt cosa 
z{a, t) = Vt sina . (4.54) 



From Eqs. ( 4.54|) , one obtains: 



pD2{a, r) = —V cosa{a? sin^ a + V'^t'^) = —V cosa\J {r'^ — a?Y + Aa?z'^ . (4.55) 

Inserting this into Eq. ( [4.1| ) and using dVt = 2tt cosa da, we have: 



V dndto y^(r2 - a2)2 + Aa^z^ 

where j^^ is the rate at which the dark matter particles fall in per unit solid angle. 
The velocity fields can also be derived. One finds: 



Vp 



— — — = ~V sina = =F — \ - ( ct^ — r^ + a/ (r^ — a^)^ + Aa'^z'^ ] 
or a V 2 V ^ / 

—1 dp'^ V'^T cos^a 



2p dr 



V 1 



^sign{z) — — \ - [r"^ — a"^ + \J(r^ — a"^)^ + Aa^z^] [r +a — \J {r"^ — a^Y + Aa^ z"^ 



where the =F signs are for the down and up flows. In the galactic plane {z = 0) we have: 

for r < a 

for r > a 
for r < a 



^ = Tvf^ 


r2 
"a2 


= 




P = 





/ a2 
= ±V\ 1 for r > a 

r 
Vu, = V- for r < a 

a 



V- for r > a . (4.58) 



r 
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V. GRAVITATIONAL EFFECTS OF CAUSTICS 

The gravitational force per unit mass caused by the particles in a zero- velocity dispersion 
flow x{d, t) is: 

d{x',t) 



Fix, t) = Gm f Sx' , y ':;, (f - X 
J \ x' — X r 



= Gm / d'a , , , « I V- ,\ - 13 ^-1 

J aaiaa2das \ x[a,t) — x l"^ 

where m is the mass of each particle. We are particularly interested in the effect of caustic 
rings on galactic rotation curves. Since caustic rings migrate only on cosmological time scales 
0, which are much longer than gas dynamic time scales, it is reasonable to expect the gas 
in the galactic plane to have relaxed to orbits consistent with the dark matter distribution 
in the caustics. 

If the gas or other material in the galactic plane at radius r moves on a circular orbit 
with velocity v{r), then 



v^(r) 



F{r) (5.2) 



is the inward gravitational force per unit mass at radius r. Let us assume that a circular 
ring caustic lies in the galactic plane at radius a. In the spirit of perturbation theory, let 
v{r) = Vrot + Vi{r) and F(r) = Fo(r) + -Fi(r), where Fi(r) is the inward force per unit mass 
due to the caustic and fi(r) is the perturbation the caustic causes in the rotation curve. If 
the matter in the caustic were smoothly distributed, the rotation curve would be flat with 

2 

value Vrot- At zeroth order, -Fo(r) = ^^^. At first order, 

2 

Fi(r) = - Vrot vi{r) . (5.3) 

r 

Let us assume that p,q,\ r — a \<< a where p and q are the transverse dimensions of the 
caustic ring. In that limit we may, when calculating -Fi(r), neglect the curvature of the ring 
and pretend that it is a straight tube. After integrating over the length of the tube, we have 

Fi(r) ~ 2Gm f dp f dz . "^^^'^^ , {r - p) ■ (5.4) 

J J [r — pY + z^ 

Using Eq. ( |4.1| ), changing variables (p, -z) — > {a, to), neglecting the (a, to) dependence of 

l^r ~ 2-Kcos(a] da^ °^^^ ^^^ ^^^^ °^ ^^^ caustic, and approximating p(a,to) by a in the 
Jacobian factor, we obtain 

^,, 2G dM /■ , , r-p(a,to) ,^ ^-. 

a dildiQ J (r — p(a, tojj + 2; (a, to) 

Ref. 1^ indicates how to extract the prefactor ^^^^^ from the self-similar infall model. Here 
we focus on the profile of fi(r) implied by the structure of caustic rings. The functions 
p(a,to) and z(a,to) are given by Eqs. ( [4.5| ) and ([4.7|) respectively (replace r by to)- One 
finds: 
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AnG dM r — a 



where ( = su/lP' and 



T(r Y\- 1 /^4 .rr. X-{T-lf + CA^ 

/(C, X) = -jdAdT ^x-iT-iy + CA^ + AA-T- ' ^'-^^ 

Fig. 8 shows I{C,^X) as a function of X for ^ = 1.0. For X < (r < a) and X > 1 (r > 
a + p), I{C,X) is constant. For < X < 1 (inside the tricusp), I{(,X) rises by an amount 
AI[(). A/ = 1 for C = 1- Fig- 9 shows A/ as a function of (. In the hmit where the tricusp 
collapses to a point {p ^ 0), there is a discontinuity in fi(r) at the caustic ring radius: 

AttG dM ^ ^_, 
A^i - 7 -^TTr^TT A/(C) • (5.8) 

VrotO d\l dlQ 

The profile shown in Fig. 8 should be added to a descending rotation curve so that the total 
rotation curve, with the effect of the caustic ring included, remains fiat on average. Fig. 10 
gives a qualitative description. 



VI. CONCLUSIONS 

We discussed the appearance of caustics in the flow of collisionless particles with negligible 
velocity dispersion. Caustics are locations in physical space where the density diverges in the 
limit of zero velocity dispersion. This happens wherever the 3D sheet on which the particles 
lie in 6D phase-space folds back. The generic caustic is a surface at the boundary between 
two regions in physical space, one of which has n flows and the other n + 2 flows. The density 
diverges as -^, where a is the distance to the surface, on the side with n + 2 flows. We 
discussed line caustics as well. These are a degenerate case where the density diverges as 
1/a where a is the distance to the line. We divided line caustics into two types: "attached" 
{Im'~f± = 0) and "isolated" {Im'-f± 7^ 0). In the latter the density is finite everywhere except 
on the line, whereas in the former the line is at the intersection of two surface caustics. 

We discussed the fall of collisionless darlc matter with negligible velocity dispersion in 
and out of a galaxy. Two types of caustic form: outer caustics which are spherical surfaces 
surrounding the galaxy and inner caustics which are rings. The caustic rings are located 
near where the particles with the most angular momentum are at their closest approach to 
the galactic center. The surface of the ring is a closed tube whose cross-section is a closed 
line which has three cusps one of which points away from the galactic center. The tube is 
the location of a generic surface caustic. Inside the tube there are four flows whereas outside 
there are two flows. In the limit where the transverse dimensions of the tube vanish the ring 
is a closed line caustic of the isolated variety. 

We analyzed the ring caustic in detail in the case of axial symmetry about the i-axis, 
of reflection symmetry z — > —z, and where the transverse dimensions, p and q, of the tube 
are much smaller than the ring radius a. The caustic is then described by 5 parameters: a, 
u, s, to and b. The precise shape of a transverse section of the ring in this limit is shown 
in Fig. 6. The 5 parameters were determined in terms of the initial velocity distribution of 
the infalling dark matter particles assuming the gravitational potential of the galaxy to be 
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spherically symmetric. Also, a qualitative description was given of the flow of dark matter 
particles on length scales of order a. 

Finally we discussed the gravitational effects of caustic rings, in particular the perturba- 
tion in a galactic rotation curve caused by a ring lying in the galactic plane. Figs. 8, 9 and 
10 describe the shape and size of the bump implied by the tricusp structure of the caustic 
ring. 
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TABLES 
TABLE I. Values of (5 as a function of j and u. 



V -^ 
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-0.1 


0.0 


+0.1 


+0.2 


ji 
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0.28 


0.29 
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0.30 


0.30 


0.10 


0.34 


0.36 


0.37 


0.38 


0.38 


0.15 


0.38 


0.40 


0.42 


0.43 


0.44 


0.20 


0.41 


0.44 


0.46 


0.48 


0.49 


0.25 


0.43 


0.46 


0.49 


0.52 


0.54 


0.30 


0.45 


0.48 


0.52 


0.55 


0.58 


0.35 


0.46 


0.50 


0.54 


0.58 


0.62 


0.40 


0.46 


0.51 


0.56 


0.61 


0.65 


0.45 


0.46 


0.52 


0.58 


0.63 


0.68 


0.50 


0.46 


0.53 


0.59 


0.66 


0.72 
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FIGURES 




FIG. 1. The wiggly line is the intersection of the {z,z) plane with the 3D sheet on which the 
coUisionless dark matter particles lie in phase-space. The thickness of the line is the primordial 
velocity dispersion. The amplitude of the wiggles in the i direction is the velocity dispersion 
associated with density perturbations. Where an overdensity grows in the non-linear regime, the 
line winds up in clockwise fashion. One such overdensity is shown. 
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FIG. 2. A generic surface caustic in phase-space (a) and in physical space (b). The two di- 
mensions {x and y) into which the caustic extends as a surface are not shown. The physical space 
density d diverges at those locations (zi and 22) where the phase-space sheet folds back. 
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a) 



b) 





c) 



d) 




e) 



f) 





FIG. 3. Infall of a turnaround sphere. The closed hnes are at the intersections of the sphere 
with the plane of the figure at six successive times. The sphere has net angular momentum about 
the vertical axis. It crosses itself between frames b) and c). After frame e) the sphere has completed 
the process of turning itself inside out. The cusps in frames d) and e) are at the intersection of a 
caustic ring with the plane of the figure. 
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a) 




C 



b) 




X 



FIG. 4. An infall sphere near where and when it completes the process of turning itself inside 
out. The lines are at the intersections of the sphere with the plane of the figure at two successive 
times, corresponding to frames c and d in Fig. 3. a labels points along the line. ^ = at points 
A and B. -^ = at point C. Points A, B and C move to point E, which is thus the location of a 



da 

caustic since D 



there. 
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FIG. 5. Qualitative description of the flow near a caustic ring. The solid lines are at the intersec- 
tion of five simultaneous infall spheres corresponding to different initial times to5 < ^04 < ••• < ^oi- 
The five numbered points are at x{a, tQk),k = 1...5, for some value of a. Points E and F are defined 
in the text. The closed dashed line is at the intersection of the caustic tube with the plane of the 
figure. There are four flows inside the dashed line whereas outside there are two. The galactic 
center is to the left of the figure. 
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p-a 



FIG. 6. Cross-section of a caustic ring in the case of axial and reflection symmetry, and 
p,q << a. 
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FIG. 7. Trajectory of a dark matter particle falling onto a galaxy, vq is the velocity at the 
initial time to- ''m is the distance of closest approach to the galactic center. 
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FIG. 8. Plot of the function /(C,^) for C = 1-0. 
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FIG. 9. Plot of the AI(C) as a function of C- 




FIG. 10. Qualitative description of the effect of a caustic ring upon a galactic rotation curve if 
the caustic ring lies in the galactic plane. The horizontal dotted line represents the flat rotation 
curve if the matter in the caustic ring is smoothly distributed. The descending dashed line is the 
rotation curve after the matter in the caustic ring has been removed. The solid line is the rotation 
curve in the presence of the caustic ring. It is the sum of the dashed curve and the profile of Fig. 
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